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I. INTRODUCTION 

Double beta decay (DBD) is a rare weak interaction process in which two identical nucleons inside 
the nucleus undergo decay with or without emission of neutrinos. The two neutrino double beta 
decay {2v(3~ (3~) which was first predicted long back by Meyer [l[ is fully consistent with standard 
model and has been observed experimentally in more than 10 nuclei. The neutrinoless double 
beta decay (fV/3 _ /3~) which involves emission of two electrons and no neutrinos, has not been 
observed experimentally and it violates lepton number conservation. The claim for observation of 
Qvf3~ (3~ decay of 76 Ge by Heidelberg-Moscow group Q is controversial and yet to be confirmed by 
other ongoing experiments. This process is one of the best probes for studying physics beyond the 
standard model. To extract the mass of neutrino via 0vj3~ (3~ decay, it is necessary to have good 
nuclear structure models for calculating the nuclear transition matrix elements (NTME) involved. 
Large number of theoretical studies for various candidate nuclei, for 2v and Qv(i~ (3~ decay, using 
many nuclear models have been carried out so that the calculated NTME can be established to be 
reliable; see for example [!, 0]. It is important to add that all models will not work for all DBD 
nuclei due to various reasons. For example for A » 64 nuclei full shell model in ( 2 }?3/2, 1 /s/2, 2 Pi/2, 
1 <7g/ 2 ) space is still not feasible. 

In contrast to the 2v j3~ j3~ decay, the positron decay modes, i.e. 2v /3 + /3 + //3 + EC/ECEC decay 
modes (hereafter, all these three combined is called 2v e + DBD) are not yet observed experimentally 
and hence there are not many theoretical studies of the NTME involved in 2v e + DBD. However, 
in the last few years serious attempts are made to measure half-lives for 2v e + DBD modes in the 
upper (pfgg/2) shell nuclei 78 Kr 64 Zn Q and 74 Se Q (in the past, attempts are also made for 
106 ' 108 Cd [|[ and 130,132 Ba [1| nuclei). Prompted by this experimental interest, recently [l(| we have 
carried out calculations for 78 Kr using the so called deformed shell model (DSM) by employing a 
modified Kuo interaction in ( 2 j>3/2, 1 /s/2: 2 Pi/2, 1 99/2) space. It is seen that the predictions of DSM 
for 2v e + DBD half-lives are close to those of QRPA and PHFB models. Extending the study in 
[Io| further, we have carried out DSM calculations for 2v e + DBD half-lives for 74 Se nucleus and the 
results are reported in this brief report. We did not consider 64 Zn as spherical shell model is well 
suited [ll[ for the three nuclei 64 Zn, 64 Cu and 64 Ni due to the fact that they are not well deformed 
with proton numbers close to the N=28 closed core. 

Over the years, we have been using with success the deformed Shell Model (DSM) based on Hartree- 
Fock states to study the spectroscopic properties, such as band structures, shapes, nature of band 
crossings, electromagnetic transition probabilities and so on, for medium heavy nuclei [13. Il3l. [lij . 
More recently this model is applied to N=Z and N=Z+1 nuclei by including isospin projection 
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[l5|,[l6[. The spectroscopic properties especially electromagnetic transitions like B(E2) and B(M1) 
values provide a stringent test for the goodness of the nuclear wave functions generated using the 
model. It is also important to add that DSM results are being used by many groups in the discussion 
of experimental data for A ~ 64-80 nuclei [17L I n addition DSM was used in calculating transition 
matrix elements for /i — e conversion in 72 Ge [181 ] and in the analysis of data for inelastic scattering 
of electrons from /p-shell nuclei [l9(. This model has also been used for studying 2v double beta 
decay transition matrix elements for 76 Ge — > 76 Se [2(| with considerable success. More recently in 
[Io| we have applied DSM to study /3-decay half lives, GT distributions, electron capture rates and 
2v e + DBD in 78 Kr. All these confirm that DSM generates good nuclear wave functions for nuclei 
in the mass region A ~ 64-80. Now, we will first discuss the DSM formalism briefly and then the 
results for 74 Se are described. 



II. DSM FORMALISM 



Half-life for the 2v e + DBD decay modes for the 0j — » Jp transitions, with Jp for the daughter 
nucleus being Jp = 0+ or 2+ , is given by 
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Ti/ 2 (k,J F )\ = G 2v (k,J F ) \M 2u (J F )\ 2 (1) 

where k denotes the modes f3 + /3 + , /3 + EC and ECEC. As, besides the 0^~ — > 0^ transition, the ECEC 
mode for 0* — > 2^~ is also of experimental interest, we are considering both Jp = 0^ and 2^ in Eq. 
(1). The integrated kinematical factors G 2v {k, Jf) are independent of nuclear structure (except for 
the dependence on the excitation ener gy E p of the Jf state of the daughter nucleus) and they can 
be calculated with good accuracy [2ll |22[ l23l. |24|. Further, the nuclear transition matrix elements 
(NTME) M 2u are nuclear model dependent and they are given by, 

VjJ+T [Eo + E N -E T } JF+1 

where |0j") , \ Jp) an d |ljv) are the initial, final and virtual intermediate states respectively and 
Ejy(Ei) is the energy of intermediate (initial) nucleus. Note that Eq = | (Ej — Ep) = \Wq. 
Here, Wq is the total energy released for different 2i/ e + DBD modes. For 0^ — > 0^ transitions 
Wo({3 + P + ) = g /3+/3+ +2m e , W (P + EC) = Q 0+EC + e b and W (ECEC) = Q ECEC -2m e + e bl + e b2 . 
Note that the Q-values are given by the difference of neutral atomic masses of parent and daughter 
nuclei involved in the positron double beta decay process and e b is the binding energy of the absorbed 
atomic electron. For the Of -> 2+ ECEC transition, denoted by ECEC*, we have W (ECEC*) = 
^ECEC — ~ 2m e + e b i + e b2 where AE is the excitation energy of the 2^ state. We have 
employed DSM to calculate the reduced matrix element appearing in Eq. ([2]). 

In DSM, for a given nucleus, starting with a model space consisting of a given set of single particle 
orbitals and effective two-body Hamiltonian, the lowest prolate and oblate intrinsic states are ob- 
tained by solving the Hartree-Fock (HF) single particle equation self-consistently. Excited intrinsic 
configurations are obtained by making particle-hole excitations over the lowest intrinsic state. These 
intrinsic states will not have good angular momentum and good angular momentum states are ob- 
tained by angular momentum projection from these intrinsic states. In general the projected states 
with same J but coming from different intrinsic states will not be orthogonal to each other. Hence 
they are orthonormalized and then band mixing calculations are performed. DSM is well established 
to be a successful model for transitional nuclei when sufficiently large number of intrinsic states are 
included in the band mixing calculations; see [l(| and references therein. Performing DSM calcu- 
lations for the parent, daughter and the intermediate odd-odd nucleus (here we need only the 1+ 
states) and then using the DSM wavefunctions, the crr~ matrix elements in Eq. (2) are calculated. 
For further details see Now we will discuss the results for 74 Se. 
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III. RESULTS AND DISCUSSION 



In our calculations of 74 Se 2v e + DBD half-lives, for the structure of the nuclei 74 Se, 74 As and 74 Ge 
we have used a modified Kuo effective interaction [25| in the ( 2 P3/2, 1 /s/2i 2 Pi/2j S9/2) space with 
56 Ni as the inert core. The single particle energies of these orbitals are taken as 0.0, 0.78, 1.08 and 
4.5 MeV respectively. DSM with modified Kuo effective interaction has been quite successfully used 
by us in describing many important features of nuclei in A ~ 60-80 region. In particular, shape 
coexistence in spectra, observed B{E2) values, band crossings and so on in 70 ' 72 ' 74 Se isotopes are 
well described by DSM [26]. Therefore, just as 78 Kr studied using DSM in [Hj, for 74 Se 2v e+DBD 
decay DSM is expected to be good. We have also verified that 74 Ge spectroscopic properties are well 
described by DSM. For 2v e + DBD half-lives calculations, we have first performed axially symmetric 
HF calculations and obtained the lowest prolate HF intrinsic states. The lowest HF single particle 
spectra for 74 Se, 74 Ge and 74 As nuclei are shown in Figs. 1, 2 and 3 respectively. Only prolate 
intrinsic states are considered in these calculations and the oblate intrinsic states are ignored just 
as in the previous 78 Kr analysis [l(| u sing DSM. The reason for neglecting the oblate states has 
been discussed in an earlier publication [27J ■ For these three nuclei, we found that the spectroscopic 
results obtained with only oblate states compare poorly with experiment and hence we did not 
include oblate states in the final calculation. We have also seen in the band mixing calculations 
that oblate states do not mix with prolate states significantly and hence they are not expected to 
affect our final results. By particle-hole excitations from the lowest intrinsic states shown in Figs. 
1-3, excited configurations are generated. For 74 Se ground state + , 10 intrinsic states with K = + 
are used for band mixing. Similarly 24 configuration with K = + for 7 Ge and 65 configurations 
with K = 1 + for 74 As are employed. We have verified that these configurations are sufficient to 
provide adequate description of 2v e + DBD. Further increase in the number of configurations does 
not change the results significantly. 

Using the wavefunctions generated by DSM, 2v /3 + EC and ECEC half-lives for 74 Se — ► 74 Ge 
transitions are calculated and the results are shown in Table I. The integrated kincmatical factors 
G2v (k, Jf) have been calculated following the prescription given by Doi and Kotani[23j]. The limits 
for /3 + EC processes in 74 Se were determined only recently, in the SuperNEMO project. Measure- 
ments of Se sample consisting of natural selenium powder using a 400 cm 3 HPGe detector resulted 
in the first Ti/ 2 limits to be > 10 18 - 10 19 yr [7] for 2v /3+EC and ECEC. The DSM results in Table 
I are the first theoretical estimates for the half-lives for positron double decay modes of 74 Se and 
there does not exist any other model calculations. Let us recall here the statement in 0]: "It is 
necessary to stress that 74 Se has never been investigated before and all results here are obtained for 
the first time. Neither has this isotope been investigated theoretically; thus there are no predictions 
with which to compare. Nevertheless, we will try to estimate the significance of the obtained results 
and the possibility to increase the sensitivity of this type of experiments in the future." 



IV. CONCLUSIONS 

In this brief report, by extending our recent results for 78 Kr [io| . we have presented results for 
positron double beta decay half lives for 74 Se. They are obtained using the DSM model with a 
modified Kuo interaction in ( 2 £>3/2, J5/2, 2 Pi/2, 1 99/2) space. As spectroscopic properties of Se 
isotopes are well described by DSM, the half-lives calculated for 2v e+DBD modes of 74 Se, given in 
Table I can be taken as reliable predictions. The calculated half-life for the ECEC mode is ~ 10 26 yr 
and it may be possible to observe this in the near future with improved sensitivity of experiments. 
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TABLE I: Experimental limit on half-lives T^ 2 along with theoretical estimates in Deformed Shell Model and corresponding phase space factor G2k for 
possible decay modes for 74 Se ^ 74 Ge. The Q- values are taken from |28l ] and Abundance (P) values are from 29]. The range (a-b) given in parenthesis 
for the theoretical estimate of the half-life is given for gA/gv ~ 1.261 and 1 respectively. 

1/2 

Transition Q-value P Decay G21/ T 2v (in yrs) 



(in keV) (in %) Mode (in yr" 1 ) Expt. limit Theory 



74 Se ^ 74 Ge 1209.7±0.6 0.89 /3+EC 2.05 x 1CT 29 > 1.9 x 10 18 [7] (14.99 - 37.9) x 10 30 

ECEC 2.63 x 10" 24 ( 7.56 - 19.12) x 10 25 

27 ^ r, r, ., m 18r^ /-, r rr on oo\ i n30 



ECEC* 3.06 x 10~ 2 ' > 7.7 x 10 ls [7] (15.55 - 39.32) x 10 



represents g.s. to 2+ state transition (AE = 595.8keV). 
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FIG. 1: HF single particle spectrum for Se. In the figure circles represent protons and crosses represent 
neutrons. The Hartree-Fock energy (E) in MeV, mass quadrupole moment (Q) in units of the square of the 
oscillator length parameter and the total K quantum number of the lowest intrinsic state are given in the 
figure. 
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FIG. 2: Same as Fig. 1 but for 74 As. 
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FIG. 3: Same as Fig. 1 but for 74 Ge. 



